TRANSPORT PROPERTIES AND NEUTRINO 
EMISSIVITY OF DENSE NEUTRON STAR MATTER 
O^ : WITH LOCALIZED PROTONS 

o\ 

D.A. BAIKO 

^ ■ A.F. loffe Physical Technical Institute, 194021, Politekhnicheskaya 26, 

St. Petersburg, Russia 
N. Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, 

00-716 Warszawa, Poland 
e-mail: baikoOastro . iof f e . rssi . ru 



OO 



> ■ AND 

^ ; p. HAENSEL 

m ' 

^^ ' N. Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, 
^^ , 00-716 Warszawa, Poland 












e-mail: haensel@camk.edu.pl 



(Received February 1, 2008) 
As pointed out by Kutschera and Wojcik, very low concentration of 

P \ protons combined with a specific density dependence of effective neutron- 

proton interaction could lead to a localization of "proton impurities" in 

^ ' neutron medium at densities exceeding four times normal nuclear matter 

density. We study consequences of the localization of protons for transport 
processes in dense neutron star cores, assuming random distribution of pro- 
ton impurities. Kinetic equations, relevant for the transport of charge, heat 

«2j ■ and momentum, are solved using variational method. Localization of pro- 

C^ ' tons removes a T~^ factor from the transport coefficients, which leads, at 

lower temperatures, to a strong decrease of thermal conductivity, electrical 
conductivity and shear viscosity of neutron star matter, as compared to 
the standard case, where protons form a Fermi liquid. Due to the local- 
ization of protons a number of conventional neutrino emission processes 
(including modified URCA process) become inoperative in neutron star 
cores. On the other hand, the energy loss rate from neutrino-antineutrino 
pair bremsstrahlung due to electron and neutron scattering off (localized) 
protons, will have a specific T^ dependence, which could modify the cool- 
ing of the neutron star core, as compared to the standard case. Possible 
astrophysical implications of the localization of protons for neutron star 
evolution and dynamics are discussed. 

PACS numbers: 97.6G.Jd, 26.60.-Hc 
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1. Introduction 



The composition of neutron star matter above three times normal nu- 
clear density (p > 3po, where po = 2.8 x 10^^ g cm~^) is largely unknown. 
Some many-body calculations suggest, that above Spo the baryon compo- 
nent of matter consists nearly exclusively of neutrons, with a small admix- 
ture (a few percent) of protons (see, e.g. models UV14+TNI, AV14+UVII 
of |l|); such a composition would be similar to that at p ~ po- A small 
admixture of protons in high-density neutron matter could behave quite 
differently than at densities p c^ pQ. As shown by Kutschera and Wojcik 
[Q], coupling of "proton impurities" to the density waves in neutron matter 
could lead, above some critical density, to localization of protons in the po- 
tential wells associated with the neutron density inhomogeneities; at smaller 
densities the coupling results in a gradual increase of the proton effective 
mass P]. We will show in the present paper, that the localization of protons 
would change dramatically transport properties of neutron star matter. It 
would also lead to the "switching on" of two new neutrino emission pro- 
cesses, which would influence cooling of a neutron star. 

In order to visualize possible effect of the localization of protons, let us 
consider the standard case, when nucleons form degenerate normal Fermi 
liquids. The transport coefficients of neutron star matter exhibit then char- 
acteristic temperature dependences, which result from the combined effect 
of the Pauli principle, and the energy and momentum conservation in the 
scattering processes y]. The leading terms in the low-temperature expan- 
sions of electrical conductivity, a, and shear viscosity, ry, are proportional 
to T~^, while the low-temperature behaviour of thermal conductivity, k, is 
given by K ex T~^. 

The conventional neutrino emission processes in the standard npe mat- 
ter include the modified URCA process and the neutrino-antineutrino pair 
bremsstrahlung from nn, np, and pp collisions. In the absence of nucleon 
superfluidity the energy emission rates for all these processes vary with 
temperature as T^ (e.g. 0|). 

Localized protons will act as scattering centers for neutrons and elec- 
trons. In what follows, we assume that the localized protons do not exhibit 
a long-range crystalline order (this point is discussed in Section 2). In such 
a case, an elementary consideration indicates, that the localization would 
give o"i°'=-P-^ ^loc.p. ^ j,o ^ ^^^ ^loc.p. ^ rp.^ Qj^g may thus expect a strong 

effect of the proton localization on the transport coefficients at lower tem- 
peratures, where scattering off protons dominates. As will be shown in this 
paper, the localization of protons produces a drastic decrease of the trans- 
port coefficients of neutron star cores, as compared to the standard case, 
when protons form a Fermi liquid: a^°'^'^'/a, fj^°'^-P-/rj ~ 10~^ — 10~^ Tg , 
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and «; P'/k ~ 10^^ Tg , for densities around Apo and for a proton fraction 
about 1 % (Tg = T/10^ K). 

Furthermore, the emission of neutrino-antineutrino pairs in the nn col- 
hsions is the only process of the above mentioned, that survives the focahza- 
tion of protons. Instead, there appear the neutrino-antineutrino pair emis- 
sion from the scattering of neutrons and electrons off localized protons. In 
both cases the rate of this process reproduces the temperature dependence 
of the direct URCA process (T^), which, if operative, accelerates drastically 
the cooling of a neutron star. As will be seen in Section 6, when neutrons are 
not superfluid the neutrino emissivity due to np collisions is approximately 
3.5 orders of magnitude larger than that due to ep collisions. Despite the 
same temperature dependence, the neutrino emission due to neutron scat- 
tering off localized protons is much less efficient than the direct URCA pro- 
cess. However, the ratio of the emissivity due to the np bremsstrahlung to 
the emissivity due to the modified URCA process (the most important one 
among the standard processes in neutron star cores) could be quite large in 
the temperature range of interest: Qg" m^ '^ /QmiiRCA ~ 2 • 10^ Tg~ , for the 
same density and proton fraction. This implies that the proton localization 
could lead to an intermediate regime of a neutron star cooling: more rapid 
than the standard cooling, provided by modified URCA process, and less 
fast than the accelerated cooling due to direct URCA process. 

The paper is organized as follows. The physical conditions in neutron 
star matter with localized protons are discussed in detail in Section 2. In 
particular, we emphasize the similarity between the behaviour of a proton in 
neutron matter and a polaron behaviour of slow electrons in solids; we also 
present some arguments against a crystalline ordering of localized protons, 
and discuss the importance of relativistic effects in the neutron component 
of matter. Kinetic equations relevant for the transport of charge, heat 
and momentum, as well as variational solutions to them, and analytical 
expressions for the transport coefficients, are derived in Section 3. Angular 
averages of scattering probabilities, appearing in the expressions for the 
transport coefficients, are calculated in Section 4. In Section 5, we improve 
the variational solutions to make them asymptotically exact in the high- 
temperature regime and present our results for the transport coefficients in 
the form suitable for practical applications. The neutrino energy emission 
rates from the ep and np bremsstrahlung are estimated in Section 6. Finally, 
in Section 7 we discuss some astrophysical implications of these results. 

Throughout the paper we will mostly use the units h = c = ks = ^ and 
turn to the normal units whenever presenting the final results. 
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Consider npe matter at super nuclear densities. The typical kinetic 
energy of a proton can be estimated from the uncertainty principle and 
reads 



n^nj^ .-(jn_ x„ \2/3 



T ^ 1112^ = 0.7 -- • ^^ \-^\ MeV, (1) 

^ 2m; V4no 0.01/ \ml j ' ^ ' 

where n is a mean nucleon density, Xp is a proton fraction (equal to an elec- 
tron fraction Xe), uq is the normal nuclear number density, no = 0.16 fm^^, 
and m* is the proton effective mass resulting from the two-body nucleon- 
nucleon (NN) interactions. 

A proton in neutron matter also has an effective potential energy T4fj, 
the value of which depends on the neutron density and ranges, according 
to different parametrization, from 55 to 75 MeV for n„ ~ 4no (fig. 1 of ||2|). 
Taking into account the possible inhomogeneity of the neutron sea we can 
write the Hamiltonian of a proton in the form 

where the last term describes the coupling of the proton to the density 
waves in neutron matter. The quantity o"(n„) = n„5T4fj/(9n„ plays role 
of the coupling strength and is of the order 550 - 600 MeV at n„ ~ 4no. 
On comparison with the estimate of the kinetic energy, Eq. (ffl) , it becomes 
clear, that the coupling to neutron density waves might not be neglected, 
if one aims at a realistic description of a proton behaviour in high-density 
neutron star matter with low proton fraction Q. 

The calculations presented in 0] indicate, that at weak and interme- 
diate couplings the proton, interacting with neutron density waves (which 
will be referred to as phonons, while, basically, these are acoustical modes 
of neutron matter), acquires an additional effective mass, which increases 
gradually with increasing neutron density. The situation reminds the so- 
called "large" polaron in solids, where a slow electron, moving through a 
crystal, is dressed into a cloud of virtual phonons, and, consequently, has 
an effective mass exceeding its "bare" band effective mass. If the coupling 
strength increases further, the polarization of the ion lattice by the electron, 
could get very strong. It can eventually be trapped by a local deformation 
of the lattice, induced by the electron itself ^. This latter situation cor- 
responds to a "small" polaron. If the temperature is sufficiently high the 
electron can be then kicked out of the potential well trapping it. However, in 
the low temperature regime, the electron could only tunnel slowly through 
the lattice and spends most of the time near one ion. 

Let us now come back to the proton moving through the neutron back- 
ground. The background polarization, induced by the proton, could be 
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characterized by its spatial scale Rp. With increasing neutron density the 
proton-phonon coupling strength a"(re„) increases, whereas the Rp decreases. 
The conditions for proton self-trapping are roughly given by two inequal- 
ities: m*(T(n„)|(5n„|i?p/n„^ > 1, which ensures, that the potential well 

— 1/3 
formed contains a bound state for a proton, and Rp < Up , which al- 
lows each proton to produce its own potential well, and ensures no overlap 
between wave functions of different protons. 

In order to check if this "small" polaron regime occurs for a proton in 
neutron matter one have to perform a detailed calculation. In particular, 
one should compare the energy of the state with trapped proton and that 
of the state, in which the proton is not localized. The results of calculations 
of this type has been reported in ref. Q. These authors evaluated the 
energy of a Wigner-Seitz cell with homogeneous distributions of neutrons 
and a proton versus that of the cell with proton wave function localized near 
the center of the cell and a neutron distribution having minimum at the 
center. The parameters, characterizing the deviation of those distributions 
from the homogeneous ones, has been treated as variational. It was found, 
that above some critical density (ranging from 4 to 9 no depending on the 
model; from now on we will adopt the most optimistic model and assume 
the critical density to be ~ 4?7-o) there existed a domain of the parameters, 
where the localized state was energetically preferable to the uniform one. 
Therefore, above this density, and provided Rp, the rms radius of the proton 
probability distribution corresponding to the minimum energy, is sufficiently 
small, the protons most likely are in effective potential wells trapping their 
wave functions. The typical value of Rp at 4no is about 1 fm and typical 
depth of the well is ~ 100 MeV. The mean distance between protons 

-1/3 = 64 t-!^. "■■ ^"'" 



while the energy of the zero-point vibration of a proton in the well is 
To r^ A2{Rp/l fm)^ {mp/m*) MeV. Thus, we see that at densities above 
4no and for sufficiently low proton fraction, say 1 - 5 %, the wave functions 
of the protons are well localized around the neutron density minima, and 
there is no overlap between them. This, in turn, means that there is no 
need to invoke Fermi statistics to describe the proton system. The temper- 
ature, when falls below, say 10 MeV, is negligible compared to the relevant 
energy scale in the proton-trapping well, so that the protons occupy the 
lowest available energy levels, and other (degenerate) particles scatter off 
them elastically. The proton could then only tunnel through the neutron 
background with very low probability, and in what follows, we will neglect 
the possibility of such events. 
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Let us discuss in some detail the astrophysical scenario of the formation 
of a dense neutron star core, where one might expect the locahzation of 
protons. According to the standard scenario of the formation of neutron 
stars in a gravitational collapse of massive stellar cores, initial temperature 
in the central core of a newly born neutron star is T ~ 30 MeV. Moreover, 
due to neutrino trapping in dense hot plasma, initial proton fraction is very 
high, Xp ~ 0.3. Both a high value of Xp and a high value of T exclude 
possibility of the localization of protons. It may take place after the neu- 
tron star core becomes transparent to neutrinos, and the core temperature 
falls to the range of < 10 MeV. One could ask a question, whether the 
proton impurities could be ordered during the localization process due to 
some long-range interaction? We immediately find, that Coulomb energy, 
estimated as, 



Ecpp ~ — ~ 0.3 (-^ ■ -^ ] MeV , (4) 



'""^^p'^^-^KMo 0.01 y 

is negligible compared to the potential energy, felt by a localized proton in 
its effective well, or the energy, resulting from the Fermi motion of protons 
before the localization point. However, it is not excluded, that some corre- 
lation of proton sites could be induced by a long-range correlation resulting 
from (strong) nuclear interactions in the neutron sea. A detailed discussion 
of this possibility is difficult and goes beyond the scope of the present paper. 
One could recall only, that the localization occurs at rather high temper- 
ature, implying large thermal fluctuations, and it is uneasy to imagine a 
long-range correlation, originating from strong (that is short-range) forces, 
with energy scale of several MeVs. Therefore, we think that neutron star 
matter with localized protons is formed as a disordered system and remains 
locked in this disordered state during subsequent cooling. In what follows, 
we will assume that the positions of localized protons show no long-range 
crystalline order: localized protons will be treated as randomly distributed, 
spatially fixed scattering centers for neutrons and electrons. 

Transport processes in neutron star matter with localized protons are 
carried out by electrons and neutrons. For simplicity, we will not allow for 
the presence of muons. Their inclusion would lead to some modifications of 
the kinetic equation formalism (see, for details, the paper by Gnedin and 
Yakovlev H, where the thermal conductivity of pe/x component of neutron 
star core matter was studied under standard assumptions) but would not 
alter qualitatively the main results of this paper. 

The Fermi energy of electrons is given by 



^- = "^-Hi^ ■ 5H I ^"=^- <=' 
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Typical energy of Coulomb ee and ep interactions are on the order of Ecpp 
(2), that is much smaller than epe- Therefore, electrons should be treated 
as a free, uniform, ultrarelativistic, and strongly degenerate {epe ^ T) 
Fermi gas. Transport of energy, momentum and charge is carried out by 
the elementary excitations in electron gas — thermal electron quasiparti- 
cles. Electron quasiparticles are approximated by the excited single-particle 
states in a free electron gas, close to the Fermi surface [i.e. with momenta 

close to kpe = (37r^ne) ], with effective mass given by ml = epe- The pair 
interaction between electron quasiparticles is described in the formalism of 
the dielectric function, the details of which are outlined in Section 3. Here 
we mention only, that since usually we have T <C ^.lepe in the neutron stars 
cores, the effect of dynamical screening can be neglected. 

Neutrons, by contrast, form a strongly interacting Fermi system. We 
will consider the case, when neutrons are normal (non-superfluid). Trans- 
port processes involving neutrons can then be described in the spirit of 
the Landau theory of normal Fermi liquids. Namely, transport of energy 
and momentum is carried out by the neutron quasiparticles — elementary 
single-particle excitations in the vicinity of the Fermi surface. As we re- 
strict ourselves to strongly degenerate neutron matter, the gas of neutron 
quasiparticles is dilute. The Fermi momentum of neutron quasiparticles 
coincides with that of the real neutron matter, kpn = (37r^n„) . The 
neutron quasiparticle velocity at the Fermi surface is given by 



Vpn = — — = 0.56 
mt 



n , 



■'='m. 



ml 



(6) 



where m* is the neutron quasiparticle effective mass. We see, that at den- 
sities of interest (we will confine ourselves to densities below 5no), and 
with Xp <C 1, neutrons gradually become relativistic, although not very 
much; nevertheless, we will take neutron relativism into account. In the 
non-relativistic system the effective mass of neutron quasiparticles m* = 
PFn/vpn differs from bare nucleon mass niN due to many-body effects. Here 
it should also include relativistic effect. Using Lorentz invariance arguments 
one obtains, in the reference frame of neutron matter, m* /e^„ = 1 + Ff/3, 
where £pn is the neutron chemical potential, which includes neutron rest 
mass, and Ff is the Landau parameter. 

Actually, the presence of proton impurities leads to the appearance of 
neutron density inhomogeneities. These inhomogeneities will slightly mod- 
ify the mean-field, felt by an incident neutron quasiparticle (the de Broglie 
wavelength of which is generally smaller than the sizes of the inhomo- 
geneities) , inducing a continuous drift of the quasiparticle in the momentum 
space. We note, that these inhomogeneities are restricted to a very small 
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fraction of the volume {Rp/ap)'^ ~ 0.01, and are concentrated around the 
proton impurities. The latter produce sudden changes in neutron momenta, 
thus making a dominant contribution to the scattering. The contribution 
from the inhomogeneities would account for a correction to the off-proton 
scattering, which, in principle, could be included into np scattering transi- 
tion probability. However, in view of the anticipated smallness of this effect, 
we will neglect it. 

3. Kinetic equation 

Transport of charge, heat or momentum is limited by scattering of the 
elementary excitations (quasiparticles) off localized protons and by their 
mutual collisions. Among the latter we neglect those caused by weak en 
electromagnetic interaction, which enables us to study transport proper- 
ties of electrons and neutrons separately. The distribution function of elec- 
trons satisfies the standard time-independent Boltzmann equation, valid for 
ideal gases, (e.g. ref. Q), whereas the neutron distribution is governed by 
the time-independent Landau equation, which takes into account a self- 
consistent mean field depending on the quasiparticle distribution itself. 
However, a proper linearization (e.g. ref. H) reveals that the quantity of 
real significance for the transport processes (the deviation of the distribution 
function from the local equilibrium one) satisfies the linearized Boltzmann 
equation, which allows a unified treatment for electrons and neutrons. We 
have 

Vl 9r/ki = ^ii + ^ipi (7) 

eE 9ki/ki = lee + lep, (8) 

where /k^ is the distribution function in question, i = {n, e}, E is an external 
electric field, vi is the velocity of quasiparticles, which is assumed to be 
independent of coordinates, and ki is the quasiparticle wavevector. On the 
right hand sides of these equations we have integrals of particle-particle and 
particle-proton collisions which can be written as 



{IttY 



Idksdk'idk^ 5{£' - 8) 5(K' - K) Cii i Y. I^-I"' (9) 



" (27r) 



(T2(T-^(T2 



2 

ip\ 



o-pO-lO-p 



/k;(l-/kJ-/k,(l-/k;)J. (10) 

In this case, non-primed and primed variables correspond to particles before 
and after a collision, respectively. An and Aip are the transition amplitudes 
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of the a and ip scattering processes, containing all exchange contributions, 
Cii is the standard "two fermion" Pauli factor, and Up is the density of 
localized protons. In the first integral the factor 1/2 serves to avoid double 
counting of the final states in a collision event. In the second one the same 
factor accounts for the fact that the density Up already includes protons 
with both possible spin orientations. 

It is useful to cast the spin sums in the preceding equations in a more 
symmetric form: 

with analogous expression for Qip. 

If the distribution functions are known, the kinetic coefficients can be 
found by calculating the currents induced in the system by the spatial gra- 
dients or the external field. We have 



JT 
je 



y'dk/k(e-M)v = -KVT, (12) 



(2^)3 

^/dk/kev = aE, (13) 

2 /■ du 

""""^ " (2^7 "^^ ^'^^"''^ " "''^' ^^^^ 

where k is the thermal conductivity, a is the electrical conductivity, r/ is the 
shear viscosity, and u is a bulk velocity, which is assumed to be directed 
along the x-axis and to be dependent only on the y coordinate (thus sat- 
isfying div u = 0). Further, Jt, je, and Uxy are, respectively, the heat flux, 
the charge flux and the xy-component of the dissipative part of the stress 
tensor. 

In the presence of the spatial gradients we approximate the true distri- 
bution functions on the left hand side of the Boltzmann equation ^ by the 
local equilibrium (le) distributions as 






e — ii'' ^ 
1 -F exp ^ 



1 -|- exp 



r(r) 

e — ^ — ku(r)"' ~^ 



(15) 
(16) 



Simple estimates Q show that it is a plausible approximation provided the 
typical scales of the gradients of temperature or bulk velocity are large 
compared to the quasiparticle mean free path /. On the left hand side of 
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Eq. (|8|) we use the global equilibrium distribution f^{£). Corresponding 
condition of applicability reads \eEl\ <C T. 

On the right hand side we adopt the following decomposition of the true 
distribution functions in terms of the f^: 

f - fO + :LilllZl$ (u) 

J K,(T,ri — J ^ rp ^K,a,ri- \^') 

In this case ^a is an unknown function of energy and angular variables, 
which represents a nonequilibrium correction to the distribution function. 
It is clear that the nonequilibrium term should be tied to the Fermi surface. 
The dependences on angular and energy variables can be separated and 
the functions $ can be looked for in the form 

% = -r,i^^(x)vV(ku). (18) 

In this case a dimensionless function \1'q represents a dependence on energy 
[x = {e — fJ')/T], while angular dependences are uniquely determined by the 
left hand sides of Eqs. (0), (g); th is a constant, which has the meaning of 
a typical time between the collisions of identical particles. 

Omitting standard details which include the linearization of the collision 
integrals Eqs. (pi), ( p!0| ) with respect to the small nonequilibrium corrections 
( p!?!) and integration of the momentum and energy conserving delta func- 
tions, we arrive at the one-dimensional equation for the functions ^a- 

5«(xi)/?(i-/?) = ^i±^/?(i-/?)*«(xi) 

+ ^f,{l-f,)il-u;^)^aixi), (19) 

Tip 

where a runs over k, a, and t], gn{xi) = xi, and ga^rj = 1- The expression 
for the characteristic collisional time is given by 

where the angle brackets denote an angular averaging of the form 

(Q,,= / ^^ft.(,,^). (21) 

J(47r) 47r COS&/2 
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In this case, the angles 9 and (p specify a coUision of identical particles 
in the Abrikosov-Khalatnikov frame of reference. Namely, is an angle 
between momenta of first and second incident particles, and (j) is an angle 
between planes, containing two incident and two final momenta. The other 
quantities characterizing ii-type collisions are 

^^ = 77^((3ii(^,</')(l + 2cos^)), 



{Qi 

\ = ^^^-(^''^)(3sin^^ sin2 </.-!)), (22) 

and Ao- = —1. The latter relation corresponds to the fact that ee collisions 
do not limit the transport of electric charge. 

A characteristic time between collisions of a particle with localized pro- 
tons reads 

_ 2 

'^ npm*kFi{Qip)' 

where in the case of zp-collisions the angle brackets stand for an averaging 
according to 

{Qip) = - r dx sin xQipix)- (24) 

vr Jo 

X is a scattering angle of a particle colliding with a proton. Finally, the 

efficiency of i]5-collisions is described by the quantities coa defined as: 

1 



Wk = w^ = -j^^^-^ {Qip{x)cosx), 

1 ,„ , . 1 



^v = JQ-^ {Q^p{x)^ (3 cos\-l)). (25) 

The Eqs. (O) can be solved exactly if either particle-proton collisions 
dominate (tjj ^ Tip) or in the opposite case. In the former situation we have 
simply ^0.(3;) = ga{x)Tip[Tii{l — u>a)]^^- If, in contrast, the zp-collisions 
are negligible the exact solution is nontrivial and can be obtained with 
the aid of a method developed by Brooker and Sykes [0] and by H0jgaard 
Jensen et al. |10|. Finally, it is possible to find variational solutions Q 



for the functions ^^ at arbitrary value of the ratio Tu/rip. This is done 
by assuming specific dependences of these functions on x (^^ ex; x, ^g-,?? = 
const), which are consistent with symmetries of Eq. (f^), and by subsequent 
determination of the unknown coefficients on integration of the equation 
over xi within infinite limits (the thermal conduction equation must be, in 
addition, multiplied by xi). The resulting expressions are: 

-1 



*«(x) 



-7^2 (3 -A,) + — (1-t^K 

Tip 



X, (26) 



12 
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'^<T,vi^) 



— IT 



:i+A.,^) + ^(1 



w, 



<^,vJ 



(27) 



When the functions ^q, are known the calculation of the transport coeffi- 
cients becomes straightforward. We insert Eqs. (|l^), (p^, (p6[), and ( p7[ 
into Eqs. (0) - (|l|) and find 



vr^ n,- T 



Ki 



(Te 



m 



3 m* {v^ii + z^^ip) ' 



e np 



[yrjii \ ^-qip 



(28) 
(29) 
(30) 



In this case UaU are the effective frequencies of particle-particle collisions: 
mfT"^ ,^ > ,„ . . 3 3 -A;, 



207r 



(Qii) (3 - Xk) = -: I'', 



rjii 



1 + A, 



(31) 



while i^aip are those of particle-proton collisions: 



■^Kip — ^f^pl^ii^Fi \S=iip) (t WkJ — Vr)ip -, 



»? 



Wk 



Wn 



■'(jep 



^Kep- 



(32) 



The above variational expressions for the transport coefficients are exact 
in the limit when ip-collisions dominate. However, in the opposite limit the 
variational method appears to be not very accurate (most notably for the 
thermal conductivity). The exact formulae of refs. |^ |l^] differs from ours 
by some factors, depending on the parameters A. We will come back to 
those factors and correct our equations in Section 5. Now we turn to a 
calculation of the angular integrals which appear in the above equations for 
the effective frequencies. 



4. Scattering probabilities 

Jf..l. Electrons 



Let us summarize briefly the main assumptions that facilitate the an- 

The 



gular averaging for electrons in Eqs. (21) - (22), and (P 
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electron-electron scattering is described adequately by a Coulomb poten- 
tial, screened by modifications of the charge density induced in the vicinity 
of each charge. This screening, being purely electron (as the protons are 
fixed and in no way respond to a small perturbation), is represented by a 
dielectric function e(w, q) [{uj, q) - is a 4-momentum transfer in a collision 
event] and results mainly in a strong suppression of the collision probability 
when the momentum transfer q is smaller than gxF = '^\/ otj-KVpe ^Fe - the 
Thomas-Fermi wave number (here a is the fine-structure constant). When 
uj <^ qv (where v is a typical velocity of electrons) the screening is static 
and the dielectric function may be approximated by ||l2| 



eO,q =1 + ^ o 7 In 

iru^ \ 3 bu 



1 — u V? 



l + u 



In 
3 



^2 



V? 



(33) 



where u = q/lkpe^ and we took into account that electrons are ultrarela- 
tivistic. In practice, the required smallness of energy transfers for a degen- 
erate system puts limit on its temperature: T <C qvpe ~ (iTFVFe = V^Tp^, 
where Tp^ is an electron plasma temperature. In terms of the degeneracy 
parameter for the ultrarelativistic system this means that T ^ Q.lepe — 
the condition we assume to be met. 

Furthermore, one can easily verify that it is reasonable to approximate 
the right hand side of the Eq. (|3^) for any q < 2kp(., by the formula 



e(0,q) = l + %, (34) 

q^ 

valid strictly only for small momentum transfers q <C 2kpe. Even with the 
latter simplification the integrations in Eqs. (21), and (p2[) remain cumber- 



some. However, for ultrarelativistic particles the ratio y = q^Ff^kFe is a 
small number and, for practical applications, it is sufficient to retain only 
the lowest order terms of expansions of the sought -for quantities in y. In 
this way we obtain 



37^3 (,4: I 

{Qee) = ^^^^, A, = --, A^ = -l + 10y2. (35) 



The resulting expressions for partial electron-electron thermal conductivity 
and viscosity coincide with those obtained by other authors (e.g. refs. [pT|,|6[) 
in the ultrarelativistic limit. 

The treatment of the electron-proton collisions is somewhat simpler. 
First of all, here, the electron screening is always static (unless, of course, 
the temperature is large enough to excite the protons at their sites). The 
situation with proton-proton correlations is not quite certain. According 
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to the arguments given in Section 2, we assume that the proton system is 
completely disordered. Under such conditions the problem of ep scattering 
is very similar to the scattering of electrons off impurities which was studied 



e.g. by Flowers and Itoh |11]. In this case the integrations can be performed 
exactly at any degree of electron relativism and at any y, but again only 
the lowest order terms in y in Eqs. (p4|), and (ESh are needed. The exact 



expressions can be found in |11], however, the formula for viscosity obtained 
in this work is slightly inaccurate. For this reason we give here both exact 
and approximate expressions for the shear viscosity and an approximate 
formula for the thermal and electrical conductivities: 

27re'^ 
iQep) = 7:^^, u;« = l + %2(lny + l), co^ = l + 12y\lny + l.25). [36) 

The exact expression for the viscosity reads 

(Qep) (1 - ^,) = -^ [1 + 2y\l + P') + 3y'l3'] In ' + ^ 



4e I y"^ 

2 



- 2-^-3yV|, (37) 



where /? = vpe/c. 

Finally, we note that the results obtained for ep scattering, remain un- 
changed, if the protons are fully spin polarized. 



4-2. Neutrons 

Contributions of neutrons and protons to the thermal conductivity and 
viscosity for the standard non-localized model of matter in neutron star 
cores were studied by Flowers and Itoh ||l3|. In the case of nucleon con- 
tributions a great deal of uncertainty is related to the description of the 
scattering processes. It is impossible to make use of the Landau theory for 
ran scattering, since at present we do not know values of the momentum 
dependent quasiparticle amplitudes in the density range of interest. Yet 
even worse is our understanding of scattering of a neutron quasiparticle at 
the Fermi surface off a localized proton. 

Another approach, fundamentally cruder, but able to supply us with 
the desired physical input, consists of a neglect of an influence of many- 
body effects on the scattering amplitudes. This means that we base our 
consideration on the data on nucleon-nucleon (NN) scattering in vacuum. 
An obvious disadvantage of this approach is an inappropriate treatment of 
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the short-range repulsive part of the NN interaction: being to some extent 
screened in a cohision of quasiparticles, it is well sampled in a collision 
of bare nucleons. Therefore, such a method is expected to overestimate 
the role of the collisional processes in limiting the neutron transport, thus 
underestimating the transport coefficients. Although, following this second 
way, we hope that it must yield correct order-of-magnitude estimate of the 
neutron transport coefficients, we warn the reader that our results can be a 
few times less than the actual values of the quantities in question. 

Within this framework the transition probabilities for nn collisions at 
any degree of relativism are easily reconstructed from vacuum differential 
cross sections in the center-of-mass (cm) reference frame, namely: 

Qnn{0,(t>) = 2 I i2 ■ 2 n / o ACr~ '^^^''^ ' ^'''^^ ' ^^^^ 

where 

^iab = -^(l-cos0), and (/>cm = </', (39) 

rriN 

are the collision energy in the laboratory (lab) reference frame and the cm 
scattering angle. 

The situation is more problematical for the case of np collisions. Under 
the condition that neutrons are not relativistic, the transition probability 
for the scattering of a neutron off a localized proton (which can be thought 
of as a neutron scattering off an external field describing by the same np 
potential) can be again derived from the differential cross sections of np 
scattering in vacuum as 

16-K'^n'^ da. 



QnpiO) = ...2 .o (-^lab, </'cm), (40) 

where 



mj^ d^r. 



2A-2 

i^iab = -^^, and 0cm = 0. (41) 

niN 

However, with growing relativism, such a procedure ceases to be adequate, 
as the scattering off an instantaneous external potential becomes different 
from the scattering of two particles (the case studied in a laboratory). But 
for the densities considered (< 5no) the neutrons are only moderately rel- 
ativistic, and, consequently, we may adopt the above formalism with the 
natural modification of the Eq. (Hol) : 



Qnp{0) = 2 ^ u2 -:o (-^lab, 0cm). (42) 



m 



+ kl dn, 



N T ^Fn 



16 D.A. Baiko, p. Haensel 



We have calculated the required angular averages of the transition prob- 
abilities [those appearing in the Eqs. (|3l|), and (^)] using the tables of 
vacuum cross sections in the cm reference frame as functions of lab energy 
(in the range from 10 to 700 MeV) and cm scattering angle. The tables 
themselves for pp and np scattering were obtained by using the partial wave 



solution WI96 available in the SAID database |14] developed at the Vir- 
ginia Polytechnic Institute by R.A. Arndt with collaborators. To construct 
the nn cross sections from the pp ones the following procedure was used. 
First of all we subtracted from the table values the well known values of the 
Coulomb cross sections. This gave us reasonable estimates of the nn cross 
sections at larger angles but made the tables inapplicable at smaller angles 
(e.g. at energies less than 400 MeV the table values at 5° became negative, 
indicating the importance of the interference terms). It was tempting then 
to extrapolate smoothly the values of the nn cross sections at larger angles 
( > 20° — 40° depending on energy) to the domain of small angles, which 
yielded a reasonable estimate of the nn cross sections over the entire an- 
gle range. Finally, at zero energy we have used the value of 3030 mb that 
followed from the theory of nn scattering length (e.g. ref. |15|| ). 

The calculations were done for the values of neutron Fermi wave vector 
kpn from 1.1 to 2.9 fm~^ (n„ from 0.5 to 5 ng). To interpolate between 
the neighbour nodes of the tables of the cross sections the bilinear interpo- 
lation was used. The results of our calculations are fitted with the mean 
and maximum errors of the fits less than 1.5 % by the following analytical 
expressions: 

^ ^ ml jQnn) (3 - A^) ^ 0.3833 o.4 

^ _ m% jQnn) (1 + A^) _ 0.1152 3.965z„ 



''" 192^7iMmb 4 2.499 + 4^' 

g ^ m% jQnp) (1 - io^) ^ 1.833 1.430^ 

""^ Wnh^lmh zl 0.3958 + z8' 



"^AT {Qnp) (1 - ^ri) _ 0.5663 4.545z„ 
16 7r?i^lmb ~ zl ^ 1.276 + 4 



_ "■'A' \-^npl y^ ^rjj _ yJ-^^^'J ^-"^"^w (^o\ 



In all these expressions the quantity Zn is the neutron Fermi wave vector in 
units of 2.666 fm~^ (corresponding to density 4no). 



5. Practical formulae 

Let us summarize the results, derived in the previous sections, and 
present them in the form useful for practical applications. First of all. 
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we note that the variational solutions obtained could be corrected to yield 
the exact asymptotes in the limit when particle-particle collisions dominate. 
This is done customarily by multiplying the partial particle-particle varia- 
tional transport coefficients by well-known factors Cq, depending on the pa- 
rameters Aq. In the case of electrons Ck(— 1/3) = 1.3, while C^(— l-|-y^) ^ 1, 
i.e. the variational solution for the shear viscosity is exact for ultrarelativis- 
tic electrons. For neutrons the coefficients Ca are weakly varying functions 
of density. For the considered density range it is a very good approximation 
to adopt fixed values of the correction factors C^ = 1.2 and C^ = 1.05. 

Bringing together Eqs. ( p8| ) - (|3^), (|35|) , (^6|), and (|4^), and incorporat- 
ing the above factors, we may write for the thermal conductivity: 



-,23 



• pro's s rm K 



V4no 0.01/ Ti 



2.2 • 10^^ — ergs s-^cm"^K"\ 



f8 



• /o pros s cm K . 



.., = 2.5.10- (^-^j n ergss-cm-K 

«n. = 8.3.1023 f^)' {^) ;! J- ergs s-^cm-^K-i, 

" li^^-olj 5^; ^^^^ s cm K , (44) 



for the electrical conductivity: 



1/3 



., = .,, = ^2.10- (^.^) .-., ,45, 



and for the shear viscosity: 



^i ^ = %^+ V' 



.ee = 1.7.10- f^.^V'^i, 

' V4no 0.01 y r| 



cm s 



1.1.1013 ("^.^L^ ^^-^^-1^-1 



r, -15. 10- ("^V ("^V^' ^ ^ g cm-^s-^ 

2 /„ n4/3 



r/n = 2.2 . 10^3 f !M^ fZ^^ 
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f^^.JtlL] _^ gcm-is-i. (46) 



n x,f 
V4^"a0iy Srjp 

Let us remind the conditions of applicability of the above expressions. 
Both electrons and neutrons must be strongly degenerate. The neglect of 
the dynamical screening effect is based on a more stringent condition for 
electrons: T <^ O.lepe. The electron number density must ensure their 
ultrarelativism. Finally, the density of neutrons is restricted to the range 
from 0.5 to 5 hq. 



6. Neutrino losses 

In this section we study two processes, contributing to neutrino cool- 
ing, which are specific for neutron star matter with localized protons. The 
first of these processes is the neutrino-antineutrino pair bremsstrahlung 
due to electron scattering off localized protons. The second process is the 
strong-interaction analog of the first one: it is the neutrino-antineutrino 
pair bremsstrahlung accompanying scattering of neutrons off protons local- 
ized in neutron medium. The temperature dependence of both these pro- 
cesses is T^, that is the same as for the direct URCA process. Therefore, 
the localization of protons may imply non-standard (accelerated) cooling 
of a neutron star. Below, we will estimate the rate of the energy losses in 
these processes under some model assumptions. 

(a) Electrons. In the case of electrons we can use the results of ref. []T6| . 
This paper was concerned with the neutrino-pair bremsstrahlung due to 
electron-nucleus scattering in the liquid phase of the neutron star crusts. 
The authors derived the following general expression for the energy loss rate 
(Eq. (8) of 0) 



Q%rem = jjj ^ iksT)%L ergs s-1 cm-3, (47) 

5o7n c* 

where Gp = 1.436 • 10~^^ ergs cm'^ is the Fermi weak coupling constant, 
factor C^ = 1.675 takes into account the emission of t'e, t'/*, and Ut, Z is the 
nucleus charge, rij is the number density of nuclei, and the quantity L sa 1, 
interpreted as a Coulomb logarithm, is a weakly varying function of Z, T, 
and rii. The authors derived also the general formula for L [their Eq. (19)], 
which takes into account the nucleus electromagnetic formfactor, the static 
structure factor of nuclei, static electron screening, non-Born corrections 
and incorporates accurately the thermal effects. They also proposed an an- 
alytical fitting formula for L [Eq. (25)] . In our case we have obviously Z = 1 



and rii = Up. The fit of [16] does not apply here (the authors were interested 
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in the crust and considered Z > 10). Besides, in our situation the temper- 
ature is much lower than the screening momentum, the low-temperature 
case of [^. In this regime the thermal effects are not important and L 
is given by much simpler formula (21) of |jl^. In this expression we can 
neglect non-Born corrections and nucleus formfactor as well as nuclei struc- 
ture factor (because our proton system is assumed to be fully disordered), 
and approximate static dielectric electron screening function by Eq. (^4[). 
Performing one-dimensional integration for ultrarelativistic electrons, we 
obtain L = 1.755. Inserting this value into Eq. (p7D we obtain 



QBt:-'- = 6.0 • 10^« T,' (£- • ^) ergs s-^ cm-^, (48) 

where Tg = T/W^ K. 

(b) Neutrons. In this subsection we will obtain an expression, which will 
enable us to estimate the energy loss rate due to neutrino pair emission from 
neutron scattering off a localized proton. To simplify the derivation we will, 
first, regard the neutrons as fully non-relativistic, and, second, will treat 
nj?-interaction in a very crude manner, assuming it to be given by a contact 
spin-independent potential, which will be treated in the Born approxima- 
tion. The interaction strength U will further be made density-dependent 
by fitting it at a given density (or, equivalently, at a given collision energy) 
to the total vacuum np cross-section, known experimentally or theoretically 
(see below). No correction of the np scattering rate to account for the pres- 
ence of the medium (except for the exclusion principle in the initial and 
final neutron states) will be made. 

The process in question, in the formalism of Feynman diagrams, is de- 
scribed (to lowest order) by two diagrams, (A) and (B). In both cases, we 
consider an initial neutron with 4-momentum k = (e, k). In the case of di- 
agram (A), this neutron first emits a neutrino-antineutrino pair with total 
4-momentum p = (w, p) = ki + k2, where ki = (wijki) and k2 = (u;2,k2) 
are, respectively, 4-momenta of neutrino and antineutrino, then propagates 
with 4-momentum k —p and finally interacts strongly with a localized pro- 
ton, absorbing 4-momentum q = (0, q) and ending in the final state with 
4-momentum k' = (e',k'). Diagram (B) corresponds to the situation, in 
which neutron first interacts strongly with a localized proton, then propa- 
gates with 4-momentum k' + p, and, finally, emits a neutrino-antineutrino 
pair of total 4-momentum p, ending in the final state k'. 

We will use nonrelativistic formalism (for neutrons) to describe the weak- 
interaction vertex, take non-relativistic neutron Green functions, when deal- 
ing with the intermediate states, and nonrelativistic spinors for neutron 
initial and final states. Then the first order matrix element assumes the 
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form 
M = Ma+Mb = -'-^ xH^^o-9AS^^ai)x [Gik-p)+Gik'+p)] /^ (49) 

where qa ~ 1-26 is the axial renormahzation constant, greek indices take 
the values 0, 1, 2, 3, and latin indices take the values 1, 2, 3, ctj are the 
standard Pauli matrices, x &nd x' are the Pauli spinors, describing the initial 
and final neutron states. Further, /^ is the neutrino neutral current given 

by 

P = ni7^(l + 75)n_2, (50) 

where ui = u{ki) and U-2 = u{—k2) are respectively neutrino and an- 
tineutrino bispinors, and bar means Dirac conjugate. In our problem the 
denominators of the neutron propagators G cannot be zero. Hence, we 
can replace them by the vacuum propagators, in spite of the fact, that our 
process goes in the presence of neutron Fermi sea. Nonrelativistic neutron 
propagator takes therefore the form 

G(u;,k) = 1— -. (51) 

uj — e]s_ + i\j 

The neutrino pair energy loss rate (in erg cm^^ s^^) is given by 

X f{l-l')uY.\M\\ (52) 



where Up is the number density of localized protons, / = /(e) and /' = f{e') 
are the neutron Fermi-Dirac distribution functions, and the summation runs 
over final and initial neutron spin states. A factor of 3 accounts for the 
emission of fe, f/x, and v-,- pairs. 

Calculating the spin-summed squared matrix elements in the straight- 
forward manner and using the identity |17| 



'^^^^ k^k',6^^Hp -k^-k2) = l (p'a'"' + '^p'pn, (53) 

UJIUJ2 o 

we obtain the following expression for the emissivity 

QB^e^- = "f^Tn' I dk dk- dp /(I - f) u: {hW + (1 - 2gi)p^] 



(2vr) 

1 
, 2mnUJ — 2pk' — p^ 2m„u; — 2pk + p^ 



X 4m2 ( ^— - " , , ) , (54) 
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where we have replaced the integration over dq by that over dp. In this 
expression we must integrate over the domain, where e' < e and for |p| < 
Lo = e — e' . In the denominators we can safely omit the terms p'^/2m,„ 
which are negligible compared to uj, and using the fact that neutrons are 
nonrelativistic, we can expand the denominators to first order in pk'/m-n 
and pk/TTi^. Then we direct the z-axis of the spherical reference frame 
along the vector k and place the vector k' into the xz-plane. Integrations 
over dp and over angle between vectors k and k' yield the emissivity in the 
form of two-dimensional "Fermi" integral 

,n-loc.p. _ oTTGjpUrip . „ r> 2 ^ lin^ 



u^ 



Qb— ----^(1 + 2.2,1)^^ 



y'dede'|k||k'|(k2 + k'')/(i-/'; 

^^"^^^'^^(l + 2.2,i)f<V4„^ 



21(27r)« ' ''^^ \m. 



'"n 



UJ^ 



dx / dcj 7 —- — , (55) 

oo io (e^ + l)(e^-^ + l)' ^ ' 

where we have taken into account that neutrons are strongly degenerate 
and have pulled all the smooth functions of the neutron momenta out of the 
integral at the Fermi surface. The latter integral is standard and is equal 
to (27r)^/504. Substituting this value into Eq. (|55|) we obtain 

The last step is to specify the value of C/^. For the contact spin-independent 
np interaction U"^ gives differential probability of elastic neutron scattering 
off a localized proton (the same for all angles 9). Our approximation consists 
in expressing U"^ in terms of the total elastic cross section for np scattering 
in vacuum, using Eqs. (^ W% - Integrating over the solid angle we have 

4.c^.l545!,«.=..fH^V ,5T) 

In this case, T^p*'^^ is the total elastic np cross section which can be obtained 
from the same theoretical model as in Section 4.2., WI96, in the SAID 
database. We propose the following fitting formula for U'^ in the energy 
range from 100 to 700 MeV: 

U' = '"^' I "^^ 5^„ where 5^, = 26.09 + ^, (58) 
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and Zn is defined below Eq. (43). The mean error of this expression is ~ 1.2% 



and the maximum error occurring at 700 MeV is ~ 2.9%. Combining Eqs. 
( |56|) and (|58|), we obtain the following numerical result: 

X r| ergss"^cm~^. (59) 



Medium effects can be expected to modify both the magnitude and the den- 
sity dependence of C/^, as compared to that given by our simple prescription, 



Eq. (58). Therefore, expression (59) should be treated as a rough estimate 
of Qr 



Brem 



7. Astrophysical implications 

As one can see from the results obtained in the preceding sections, pres- 
ence of randomly distributed localized protons, removing the T~^ factor 
from the transport coefficients, leads at lower temperatures to a dramatic 
decrease of k, a, and rj, as compared to the standard case, where protons 
form a Fermi liquid. Let us discuss in some detail the astrophysical impli- 
cations of this effect for neutron stars. 

(a) Thermal conductivity. In our case transport of heat is dominated 
by electrons, which scatter predominantly off localized protons, so that 
^loc.p. ^ k|P'^-p- ~ k}°^'P\ Let us compare this result with the values of 
K corresponding to the standard case. The latter are expressed by the for- 
mulae derived in |p: 

-e = —^^l^-j [-) ergs cm s K (60) 

for non-super fluid p, and 

5.6-102Vne\ -1 -iT^-i 



^ , — ergs cm s K (61) 

for highly superfluid p, the decisive difference between the two regimes being 
the presence or absence of proton screening. At n = 4no and Xe = 0.01 we 
thus get 

Ke (5 ^ 10) • 10^ 



jloc.p. 2^2 



(62) 



The diffusive thermal conductivity due to neutrons is generally on the order 
of or somewhat smaller than that of electrons. (The actual ratio between 
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them depends sensitively on i) the fraction of electrons, ii) the nucleon 
effective masses, and iii) the strength of quasiparticle interactions; in the 
non-super fluid regime k„^'''' may serve as a good approximation of the neu- 
tron conductivity.) Thus we can accept the value of 10^/Tg as a reasonable 
estimate of the ratio k/k P' . 

Neutron stars are born as very hot objects which cool subsequently due 
to neutrino losses from their interior. The cooling is accompanied by thermal 
equilibration of the stellar interior with a typical time scale, determined by 
the size of the core Rcore, the specific heat per unit volume C, and the 
thermal conductivity of matter k as 

Tt.e. ~ • (63) 

K 

For estimates, in the standard non-superfluid case we can set C = Cn, the 
specific heat of degenerate neutrons, 

Cn = 2.6 • lO^o (^) (^) '^' Tg ergs cm-^K^S (64) 

and approximate k by k^ (pO). Approximating further all the effective 
masses by bare masses and assuming a core of constant density {n c::^ Un = 
4no) and composition (xg = 0.01) we obtain Tt.e. ~ 460 Tg (i?core/10 km)^ 
years. If the core is superfluid we may use for C the specific heat of electrons 

Ce-. 

Ce = 6.6 • 10^8 (J^ . J^\ ' ' j.g ergs cm-^K-i, (65) 



V4no 0.01, 

and take k from Eq. (|6l|). This gives us rt.o. ~ 10 Tg (i?core/10 km)^ years 
(the latter relation should be regarded as a lower bound, as the neutron spe- 
cific heat could decrease rather slowly depending on the type of superfiuidity 
]l8[). The temperature dependence of both those expressions is removed by 
the localization of protons. Multiplying Tt.e. by the factor k/k P' (p^) we 
obtain t^^;^' equal to 4.6-10^ {R^°^^i- /IQ km)^ years and 10^ {R^°^^i-/IQ km)^ 
years for non-superfluid and strongly superfluid cases, respectively. 

The thermal conduction becomes important when the above time scales 
are comparable or smaller than a characteristic time of thermal evolution 
of matter due to neutrino losses. Otherwise, the temperature of a given 
element of matter is determined locally. The neutrino cooling time scale 
could be estimated as CT/sy, where e^ is the total neutrino emissivity. 
The temperature dependences of e^ for various neutrino-emission processes 
are usually very strong. In view of this, the equilibration temperature for 
the case of proton localization will not be much smaller, than that for the 
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standard case; in both cases the temperatures fah to the range of 10* - 10^ 
K. However, the time required to reach the thermal equihbrium in a core, 
containing locaUzed protons, will be two orders of magnitude longer than 
that for a standard liquid core. The latter conclusion will not be changed 
by the appearance of neutron superfluidity. 

Another difference between the two cases comes from the simple idea 
that the thermal equilibration time in the localized protons case is tem- 
perature independent. While in the standard case the core below ~ 10* 
K could be treated as perfectly isothermal, it will still take ~ 10'^ — 10^ 
years to wash out any accidental temperature inhomogeneity in the core 
with localized protons. 

(b) Electrical conductivity. The electrical conductivity is relevant for the 
ohmic dissipation of internal magnetic fields in neutron stars. In the stan- 
dard case of npe matter with non-super fluid protons the charge transport 
is dominated by ultrarelativistic electrons, scattering off protons, which re- 
sults in an electrical conductivity a^p = 2.1 x 10^^ {rie/no)^''^ /Tg s~^ [19|. 



Localization of protons prevents the appearance of a proton superconductor 
(as the localization temperature is rather high, of the order of 10^'^ — 10^^ 
K, which is commonly believed to be higher than Tcp) and at n = 4no and 
Xe = 0.01 reduces the electrical conductivity of neutron star matter by a 
large factor 



aep 2 • 10^ 



loc.p. '~r2 



(66) 



Such a low value of electrical conductivity would lead to a significant decay 
of the electric currents, circulating within the core with localized protons, 
over a time scale of 

loc.p. ('r>loc.p.N2 / oloc.p. \ 2 

loc.p. ^ ^e l^corc ) ^ 3 . iqM tls^^ years. (67) 

c^ \ 10 km / 

Current analyses of the population of radio pulsars do not show any 
evidence of magnetic field decay during active lifetime of a normal radio 
pulsar ||2^. Specific lower hounds on the magnetic field decay timescale, ob- 
tained using various types of statistical analyses, range from 2 x 10^ years to 



10* years |2C]. To be consistent with this observational fact, we should as- 
sume that either the core with localized protons is small, and the bulk of the 
field is sustained by the currents, circulating in the outer part of the core; 
either the external field is separated from the internal field, and the observ- 
able bounds for the decay of the surface field put no evident constraints on 
the evolution of the core field; or that the core magnetic field is due not to 
electric currents but results from a permanent ferromagnetic magnetization 



of the matter. Actually, as demonstrated by Kutschera and Wojcik [21, 22 1, 
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ferromagnetism due to a complete spin polarization of protons and a partial 
spin polarization of neutrons (of the order of Xp) could be a generic property 
of the neutron star matter with localized protons (for a confrontation of this 
theoretical prediction with observations of radio pulsars, see [p^). 

(c) Shear viscosity. The localization of protons leads also to a strong 
decrease of the shear viscosity of neutron star matter. Assuming normal 
neutrons, we can estimate the value of r] of standard npe matter by the 
quantity r/J^^P' . The localization of protons will result in a decrease of rj by 
a factor 

ri 7-105 



^loc.p. 2| 

In contrast to the standard case yyio'^P- is temperature independent. This 
might be relevant for stability of rapidly rotating neutron stars. In the 
standard case of the npe matter, r] increases with decreasing temperature 
as T~^. Dissipative effects due to rj contribute to damping of the secu- 



lar instability driven by the gravitational radiation reaction (GRR) [p3, 25 1 
in rapidly rotating neutron stars. Detailed calculations show, that viscous 
effects of r] damp completely the GRR secular instability, if internal tem- 
perature falls below 10^ K ||2^. However, within a neutron star core with 
localized protons, the shear viscosity remains constant and close to the value 



of the shear viscosity of standard npe matter at T ~ 10"^^ K [see Eq. (38)]. 
Such a low value of the shear viscosity could not prevent the growth of the 
GRR secular instability at any internal temperature of a neutron star [p^]. 
(d) Neutrino cooling. At the earlier stages of the evolution matter 
with localized protons cools emitting neutrinos via two reactions involv- 
ing only nucleons, and several reactions, which involve also electrons. The 
two nucleonic reactions are the neutrino-antineutrino pair bremsstrahlung 
in the nn and np collisions. The first process (e.g. Q) is common for 
both standard and localized protons models of matter, and the rate of 
energy emission in this process varies with temperature as T^. The sec- 
ond nucleonic process is modified drastically by the localization of pro- 
tons: its emissivity becomes proportional to T^, reproducing the tem- 
perature dependence of the direct URCA process. Nevertheless, the np 
bremsstrahlung remains several orders of magnitude less efficient than the 
direct URCA (Qg" m^ '^ /QdURCA ~ 3 • 10"^). On the other hand, at tem- 
peratures below ~ 10^ K its emissivity exceeds that of the modified URCA 
(Qsrem'^ '^ /QmURCA ~ 2 • 10^ T^'^) process, which, in turn, is thought to 
govern the cooling in the standard (non-localized) model, when the direct 
URCA is forbidden by the momentum conservation law. Thus, one might 
expect that cooling of a neutron star with significant fraction of mass being 
in the phase with localized protons would follow some intermediate path be- 
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tween the curves describing the standard (modified URCA) and accelerated 
(direct URCA) coohng. 

As is any nucleonic process, the np bremsstrahlung is subject to strong 
suppression, if the neutron superfluidity appears. In this case, the neutrino 
pair emission from the scattering of electrons off localized protons serves as 
a dominant mechanism of cooling of a neutron star core. 

(e) Final remarks. In this subsection we will come back to the problem of 
a crystalline ordering of localized protons and comment on the values of the 
transport coefficients in this case. As discussed in Section 2, we consider this 
possibility unlikely, however, see ref. |2^. For temperatures below 1 MeV 
the phase space available for phonons (in the proton crystal) is small (as 
T^). Then the transport of energy, charge and momentum would be limited 
by ee, nn, and, in the case of the electrical conductivity, by en collisions. 
This means that all the transport coefficients would be even larger than 
those in the standard case and would reproduce the standard (Fermi liquid) 
temperature dependences. The situation would be complicated by the band 
structure of quasiparticle states, however, this would further increase the 
transport coefficients by reducing the phase space available for scattering 
quasiparticles. 
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